PL-TR-96-2105 

Environmental  Research  Papers,  No.  1192 


REPORT  ON 

PHASE  BLOCKAGE  PHENOMENOLOGY  SESSION, 
SSA  91  ST  ANNUAL  MEETING, 

2  APRIL  1996,  ST.  LOUIS,  MISSOURI 


Delaine  T.  Reiter  and  Rong-Song  Jih 


5  JUNE  1996 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


PHILLIPS  LABORATORY 
Directorate  of  Geophysics 
AIR  FORCE  MATERIEL  COMMAND 
HANSCOM  AFB,  MA  01731-3010 


19970117  078 


miG  qUiJJTY  i 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


JAM^.  LEWKOWICZ 
DireclT 

Earth  Sciences  Division 


7 


This  report  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is  releasable  to 
the  National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical  Information 
Center.  All  others  should  apply  to  the  National  Technical  Information  Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing  list,  or  if  the 
addressee  is  no  longer  employed  by  your  organization,  please  notify  PL/IM,  29  Randolph 
Road,  Hanscom  Ara,  MA  0173 1-3010.  This  will  assist  us  in  maintaining  a  current 
mailing  hst. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a  specific 
document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  re^onse,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1216  Jefferson 
Davis  Highway,  Suite  1  204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20603. 


1 .  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


REPORT  DATE 

5  June  96 


3.  REPORT  TYPE  AND  DATES  COVERED 


Report  on  Phase  Blockage  Phenomenology  Session,  SSA  91st  Annual  Meeting,  2 
April  1996,  St.  Louis,  Missouri 


6.  AUTHORIS) 

Delaine  T.  Reiter  and  Rong-Song  Jih 


1  Apr  -  30  Apr  96 


5.  FUNDING  NUMBERS 

PE  621010 
PR  7600 
TA09 
AVU08 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Phillips  Laboratory/  Earth  Sciences  Division  (GPE) 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

PL-TR-96-2105 
ERP,  No.  1192 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS{ES) 

HQ  AFTAC/TTR  Phillips  Lab  (GPE) 

1030  South  Highway  AlA  29  Randolph  Road 

Patrick  AFB,  FL  32925-3002  Hanscom  AFB,  MA  01731-3010 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12a,  DISTRIBUTION  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


12b.  DISTRIBUTION  CODE 


13,  ABSTRACT  (Max/mum  200  words) 

This  report  details  the  results  of  a  special  session  on  regional  phase  blockage  which  was  held  during  the  Seismological 
Society  of  America's  91st  Annual  Meeting  in  St.  Louis,  Missouri  on  April  2,  1996.  The  session  consisted  of  seven  orally 
presented  papers  by  members  of  the  seismic  monitoring  research  community  and  a  follow-up  meeting  that  was  intended  to 
be  a  discussion  session  about  the  scientific  problem  of  phase  blockage.  The  report  contains  the  abstracts  from  the  talks, 
background  on  regional  phase  blockage,  some  results  from  the  papers  presented  at  the  session,  and  some  conclusions  and 
recommendations  for  future  research  by  the  authors. 


14.  SUBJECT  TERMS 


Seismic  wave  propagation.  Phase  blockage.  Regional  seismic  monitoring 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


1 7.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev.  2-89)  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


CONTENTS 


SUMMARY  .  1 

BACKGROUND  INFORMATION  .  3 

SIGNIFICANT  FINDINGS  FROM  THE  SSA  MORNING  SESSION  .  5 

CONCLUSIONS .  6 

RECOMMENDATIONS .  7 

REFERENCES  .  9 

APPENDIX:  ABSTRACTS  PRESENTED  AT  PHASE  BLOCKAGE  SESSION  .  11 

DISTRIBUTION  LIST  .  19 


ill 


IV 


PREFACE 


This  project  was  conducted  under  the  auspices  of  the  Air  Force  Technical  Applications 
Center  and  the  Earth  Sciences  Division  of  the  Air  Force  Phillips  Laboratory.  The  views  and 
conclusions  contained  in  this  report  are  those  of  the  authors  and  should  not  be  interpreted  as 
representing  the  official  policies,  either  expressed  or  implied,  of  the  U.S.  Air  Force  or  the  U.S. 
Government. 


V 


j 

! 


Vi 


SUMMARY 


A  spedal  session  on  Regional  Phase  Blockage  Phenomenology  V4&s  held  on  April  2,  1996,  at  St. 
Louis,  Missouri,  during  the  SSA  91st  Annual  meeting.  The  session  was  a  result  of  several  discussions 
between  Phillips  Laboratory  (PUGPE),  AFOSR,  and  DOE’s  national  laboratories  after  the  17th  Seismic 
Research  Symposium.  The  motivation  was  to  gather  Comprehensive  Test  Ban  Treaty  (CTBT)  research¬ 
ers  to  review  the  status  and  current  direction  of  research  on  blockage-related  topics,  and  if  appropriate, 
make  some  suggestions  on  follow-up  topics.  The  morning  session  of  oral  presentations,  presided  over 
by  Dr.  Delaine  Reiter  (PL/GPE)  and  Rong-Song  Jih  (PL/GPE),  was  followed  by  a  wrap-up  discussion 
that  evening  convened  by  Dr.  Stanley  Dickinson  (AFOSR)  and  Dr.  Delaine  Reiter  (PL/GPE).  This  report 
provides  some  bad^round  on  regional  phase  blockage  and  summarizes  the  results  from  the  morning 
session  and  the  issues  that  were  addressed  at  the  discussion  session.  The  seven  abstracts  presented 
at  the  morning  session  are  induded  in  an  Appendix  to  this  report  with  titles  of  the  talks  listed  below; 

1.  Baumgaidt,  Lg  props^ation-path  barriers  in  the  Eurasian  continental  craton  -  possible  shallow  crust 

explanations. 

2.  Priestley,  Schultz,  and  Patton,  Anomabus  surface  wave  propagation  across  the  south  Caspian  basin 

and  the  blockage  of  regional  seismic  phases. 

3.  Jih,  McLaughlin.  Dainty,  and  Harkrider,  Numerical  modeling  of  Rg  blockage  and  comparison  of  Rg-to- 

Lg  conversion  mechanisms  due  to  scattering  and  incomplete  dissipation. 

4.  Xie,  Pn  amplitude  variation  aloiy  the  northern  foot  of  the  Tienshan  mountains. 

5.  Anderson  and  Cormier,  Regional  variations  in  Lg  observed  and  synthesized  at  the  CNET  and  KNET 

arrays. 

6.  McNamara,  Walter,  Schultz,  and  Goldstein,  Regional  phase  propagation  in  northern  Africa  and  the 

Mediterranean. 

7.  Sweeney,  Interpretskion  of  crustal  phase  characteristics  in  Iran  and  the  surrounding  region  determined 

from  ILPA  data. 

These  papers  will  be  referred  to  by  their  title  numbers  throughout  the  remainder  of  this  report. 
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BACKGROUND  INFORMATION 


The  most  prominent  regional  ii^ase  generated  by  shallow  events  is  the  Lg  phase  identified  by 
Press  and  Ewing  (1952).  The  Lg  phase  has  received  considerable  attention  in  recent  years  because  of 
its  potential  use  in  yield  estimation  {e.g.,  Nuttli,  1986)  and  discriminating  between  earthquakes  and 
explosions  at  regional  distances  {e.g.,  Blandford,  1981;  Pomeroy  et  al.,  1982).  Based  on  empirical 
observations,  Lg  has  been  recognized  as  a  stable  relative-yield  indicator  (Patton,  1988;  Hansen  et  al., 
1990).  However,  it  is  also  recognized  that  Lg,  like  Rg,  is  sensitive  to  changes  in  structure  along  its  path, 
which  can  have  deleterious  effects  on  any  role  as  a  discriminant  {cf.  Lynnes  and  Baumstark,  1991)  or 
magnitude  measure  (cf.  Jih  and  Lynnes,  1993),  unless  the  propagation  effects  are  accurately 
accounted  for.  Numerous  studies  have  used  the  sensitivity  of  Lg  to  structural  effects  to  map  regions  of 
anomalous  propagation  and  try  to  associate  them  with  crustal  structure.  Bias  in  Lg-based  magnitude 
measurements  were  reported  by  Gregersen  (1984)  in  Greenland.  Very  low  Lg/S„  amplitude  ratios  have 
been  observed  after  crossing  the  Tibetan  plateau  (Ruzaikin  et  al.,  1977),  the  North  Sea  grabens  (Gre¬ 
gersen,  1984;  Kennett  et  al.,  1985),  the  Caspian  Sea  or  the  Black  Sea  (Levshin  and  Berteussen,  1979; 
Kadinsky-Cade  et  al. ,  1981).  Chinn  et  al.  (1980)  (and  Lynnes  and  Baumstaik,  1991)  observed  that  the 
efficiency  of  the  Lg  propagation  is  better  for  paths  parallel  to  the  structural  trend  than  for  paths  in  the 
perpendicular  direction.  Earlier  studies  (Press  and  Ewing,  1952;  staikr  and  Ewing,  1957;  Savarensky 
and  Valdner,  1960)  have  established  that  Lg  does  not  propagate  through  crust  overlain  by  water  deeper 
than  2  km.  On  the  other  hand,  propagation  across  a  marginal  sea  of  continental  shelf  does  not  com¬ 
pletely  quench  Lg,  but  can  reduce  its  amplitude.  Baumgardt  (1991)  compared  the  crustal  cross-sections 
for  Lg  propagation,  and  found  that  the  Lg  blockage  correlates  with  thick  sediment  very  well.  Basically,  his 
observation  was  that  paths  that  do  not  cross  basins  or  for  which  sediments  do  not  vary  by  greater  than 
3  km  exhibit  little  or  no  Lg  blockage  and  scattering.  Baumgardt  (1991)  also  identified  paths  for  which  the 
surface  elevations  and  crustal  thickness  change  substantially,  and  yet  Lg  propagates  effidently.  Thus  his 
observations  suggest  that  the  near-surface  sediment-thickness  variations  seem  to  correlate  more 
strongly  with  Lg  blockage  than  do  the  crustal-thickness  variations.  Zhang  and  Lay  (1994)  used  surface 
topography  as  a  manifestation  of  the  varying  crustal  structures.  They  found  a  strong  correlation 
between  S„/Lg  ratios  for  Eurasian  explosions  and  roughness  or  mean  altitude  of  the  topography  along 
the  path,  based  on  a  meager  data  set. 

Despite  the  long-time  interest,  many  fundamental  questions  about  the  excitation  and  propagation 
of  Lg  waves  phases  remain  to  be  answered.  Numerical  modeling  of  the  Lg  waves  (and  other  regional 
phases)  would  complement  the  empirical  studies  by  providing  more  accurate  interpretations  and  better 
insight  of  the  underiying  physics.  Theoretical  studies  of  Lg  propagation  across  continental  margins  have 
been  conducted  by  Kennett  (1986),  Maupin  (1989),  Regan  and  Harkrider  (1989),  Cao  and  Muirhead 
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(1993),  and  Gibson  and  Campillo  (1994),  using  different  techniques.  A  simple  geometrical  ray  theory 
can  be  used  to  predict  the  kinematic  property  in  a  qualitative  manner,  as  Kennett  (1986)  has  illustrated, 
but  would  fail  to  explain  the  dynamic  properties  for  complex  media.  Kennett  and  his  associates  used  a 
modal  summation  to  investigate  Lg  propagation  in  stratified  and  weakly  heterogeneous  media  Mitchell 
and  Hwang  (1987)  computed  multi-mode  synthetics  for  1-D  models  with  various  thicknesses  of  low-Q 
sediments.  Regan  and  Harkrider  (1989)  used  a  hybrid  of  propagator  matrix  and  finite-element  [FE] 
methods  to  model  the  SAY-Lg-wave  propagation.  Cao  and  Muirhead  (1993)  applied  a  2-dimensional  P-SV 
finite-difference  method  to  explore  Lg  blockage  and  argue  that  a  water  column  over  the  crust  is  an 
important  factor  in  blocking  Lg  propagation.  Gibson  and  Campillo  (1994)  applied  both  the  dynamic  ray 
tracing  and  the  boundary-integral  equation  methods  to  model  Lg  blockage  in  the  west  Pyrenees  Range, 
near  the  French-Spanish  border.  They  suggest  that  the  unmod^ed  scattering  by  small-scale  features 
within  the  lower  crust  is  the  reason  for  the  observed  blockage.  Jih  (1995, 1996b)  implemented  a  pure  Lg 
wave  packet  for  Linear  Finite  Difference  (LFD)  calculations.  This  strategy  works  particularly  well  for 
models  with  additional  levels  of  complexity  due  to  free-surface  topographic  irregularities  and  anelastic 
attenuation.  Jih  (1996a)  (and  [3])  pointed  out  that  a  shallow,  strong  attenuative  layer  near  the  surface 
could  cause  Rg-to-higher  mode  coupling.  As  a  result,  some  undissipated  Rg  energy  could  propagate  as 
crustal  S  waves  like  Lg.  This  is  rather  surprising  and  seems  to  be  a  special  mechanism  only  affecting  Rg. 
The  significarx^e  of  the  LFD  modeling  research  by  Jih  and  his  partners  at  both  former  Teledyne  Geotech 
and  Phillips  Laboratory  is  that  their  distinct  code  permits  a  convenient  investigation  of  the  path  effect 
(due  to  a  specific  mechansim)  on  a  specific  phase. 
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SIGNIFICANT  FINDINGS 
FROM  THE  SSA  MORNING  SESSION 


From  the  seven  talks  from  the  morning  session,  there  are  several  observations  worth  reporting: 

[A]  With  the  exception  of  one  theoretical  study  [3],  all  the  other  talks  were  based  primarily  on  empirical 
observations.  A  few  of  the  empirically-based  talks  were  supplemented  with  minor  modeling  efforts. 
The  synthetic  techniques  used  include  the  LFD  method  (e.g.,  [2]),  ray-tradng  (e.g.,  [5]),  and  a 
hybrid  method  using  LFD  and  locked  mode  (e.g. ,  [3]). 

Several  of  the  authors’  results  indicate  that  shallow,  lateral  crustal  heterogeneity  can  explain  much 
of  the  variability  in  high-frequency  Lg  and  long  period  surface  waves  in  the  continents  ([1-3], [6]).  At 
least  two  of  the  studies  indicate  that  sedimentary  basins  combined  with  thinned  crust  are  espe¬ 
cially  efficient  at  attenuating  Lg;  such  structures  are  easily  modeled  using  2D  approaches.  One 
question  that  then  arises  is  how  geological  structures  that  are  not  sedimentary  basins  can  block  Lg 
([1]  used  the  example  of  the  Tsangpo/Indus  suture  zone),  and  how  these  effects  could  be 
modeled. 

The  principal  conclusion  reached  by  the  authors  of  [^  is  that  strong  gradients  in  Moho  topography 
combined  with  increasing  length  of  path  across  such  topography  correlates  well  with  decreased  Lg 
efficiency.  This  contrasts  with  the  results  from  papers  [1],  [2]  and  [4],  which  indicate  that  Moho 
topography  is  a  secondary  effect  compared  to  upper  crustal  heterogeneity.  One  possible  explana¬ 
tion  for  this  discrepanq^  is  that  the  modeling  technique  (ray-tradng)  may  bias  the  author’s  condu- 
sions  towards  topography  as  a  predominant  cause  of  surface  wave  blockage. 

While  the  Lg  phase  remained  the  theme  of  several  talks,  several  Pis  have  investigated  other  regional 
phases  such  as  Rg  and  Pn.  This  kind  of  study  is  important  because  blockage  of  one  phase  in  the 
crustal  and/or  mantle  waveguide  is  often  assodated  with  coupling  or  conversion  to  another  phase. 
For  instance,  in  [3]  Jih  et  al.  numerically  modeled  pure  Rg  using  2D  linear  finite  difference  (LFD) 
and  compared  a  number  of  physical  mechanisms  and  their  effects  on  waveforms.  These  mechan¬ 
isms  induded  rough  topography  and  random  heterogeneity.  They  also  looked  at  Rg->Lg  conver¬ 
sion  in  the  presence  of  shallow,  highly  attenuative  layers.  The  relative  effectiveness  to  block  Rg 
using  all  three  mechanisms  was  compared.  This  study  provided  a  dear  picture  of  the  expected 
effect  on  waveforms  using  three  important  and  viable  physical  mechanisms. 

In  [4]  Xie  used  Pn  amplitude  variation  as  an  indicator  of  crustal  heterogeneity.  He  found  that  Pn 
variation  is  much  more  unpredictable  than  other  phases.  It  is  also  more  difficult  to  tell  what  is 
causing  the  variability;  it  could  be  a  number  of  mechanisms  such  as  3-D  heterogeneity,  caustics  at 
the  surface  and  Moho,  or  near-surface  structural  variation.  Further  study  is  necessary  to  quantify 
which,  if  any,  of  these  mechanisms  can  explain  the  variability  in  Pn. 

[C]  Many  of  the  presented  talks  utilized  a  variety  of  seismic  data  sources  ranging  from  recently 
recovered  historical  ILPA  data  ([7])  to  data  recorded  at  the  newly  deployed  KNET  and  CNET 
([4,5]).  The  significance  of  utilizing  these  new  data  sources  in  regions  of  great  CTBT  concern  is 


4 


that  it  demonstrates  the  reliability  and  usefulness  of  the  new  networks.  It  is  hoped  that  further  stu¬ 
dies  using  data  from  regions  of  CTBT  interest  will  further  complete  maps  of  regional  phase  block¬ 
ages  and  explain  some  of  the  mechanisms  which  cause  blockages. 
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CONCLUSIONS 


Lg  phase  blockage  on  the  continents  appears  to  be  dosely  tied  to  upper  crustal  structure.  The 
modeling  and  empirical  studies  show  that  Lg  waves  traversing  regions  with  thick  sedimentary  basins  or 
other  strong  lateral  variations  are  effectively  blocked.  A  second  and,  to  our  minds,  less  important  cause 
of  phase  blockage  is  that  of  topography,  both  surface  and  Moho.  Other  phases,  such  as  Pn,  are  also 
sensitive  to  upper  crustal  structure  and  Moho  topography,  but  they  are  less  well  understood. 

If  the  empirical  studies  are  done  correctly  (that  is,  they  separate  out  the  effects  on  blockage  of 
source  mechanism  and  other  forms  of  attenuation  from  that  due  to  pure  path  blockage  from  topographic 
or  structural  variations),  they  can  provide  a  valuable  database  or  "map"  of  blockage  in  regions  of  CTBT 
interest.  An  empirical  approach  serves  the  research  goal  well  to  the  extent  of  confirming  the  phase 
blockage.  However,  once  the  blockage  is  observed  in  a  new  region  (or,  in  some  cases,  a  well-studied 
region),  perhaps  the  effort  should  be  directed  towards  improving  the  fundamental  understanding  of  (a) 
the  physical  mechanism  as  well  as  (b)  the  possible  implication  for  event  discrimination.  These  two 
branches  generally  reflect  the  difference  between  6.1  and  6.2  projects.  Obviously  from  a  CTBT- 
monhoring  point  of  view,  it  is  never  enough  to  "re-confirm”  a  phase  blockage  phenomenon.  At  some 
point,  the  Pi  should  move  fonvard  to  either  (a)  try  to  solve  the  sdentific  puzzle  or  (b)  find  some  correc¬ 
tion  procedure  to  account  for  the  phase  blockage.  Either  one  will  improve  our  capability  in  monitoring  a 
CTBT  in  the  regional  regime.  Very  often,  though,  these  two  goals  are  intertwined  and  difficult  to 
separate. 

Modeling  studies  done  in  conjunction  with  empirical  observations  (such  as  papers  [2]  and  [6])  can 
reveal  the  causes  of  blockage,  as  long  as  they  use  an  appropriate  modeling  technique  that  can  simulate 
the  full  wavefield.  Studies  like  [3]  are  important  to  understanding  the  physical  mechanisms  of  a  particu¬ 
lar  type  of  propagation  and  blockage,  in  this  case  Rg. 
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RECOMMENDATIONS 


In  our  opinion  the  problem  of  phase  blockage  divides  fairly  naturally  into  6.1  (basic)  and  6. 2/6.3 
(applied)  type  research  programs.  A  basic  research  program  should  be  focused  on  adding  to  our  further 
understanding  of  the  physical  mechanisms  of  Lg,  Sn  and  Pg/Pn  blockage.  This  includes  work  on  estab¬ 
lishing  the  physical  basis  for  discriminants  in  terms  of  source  mechansm  and  regional  phase  excitation. 
Without  a  firm  foothold  in  source  theory,  questions  about  transportability  will  continue  to  go  unanswered, 
particularly  in  regions  where  few  if  any  nuclear  explosions  have  been  tested  before. 

Establishing  the  physical  basis  can  be  very  difficult  because  there  are  usually  several  competing 
mechanisms,  which  can  often  be  coupled  and  poorly  understood.  A  particularly  challenging  puzzle  con¬ 
cerns  the  physical  basis  for  the  P-to-S  ratio  and  the  spectral  ratio  discriminants.  The  focus  here  is  pri¬ 
marily  on  the  generation  of  explosion  Lg  waves  and  the  role  of  spall  versus  the  role  of  Rg-to-S  conver¬ 
sions. 


Rg-to-S(Lg)  conversion  in  the  source  region  is  a  propagation  effect  involving  a  disrupted  waveguide 
or  heterogeneity  in  the  medium.  On  the  other  hand,  explosion  Lg  may  be  generated  predominantly  by 
intrinsic  source  mechanisms,  especially  at  high  frequencies.  (See  McLaughlin  et  al.,  1988;  Patton  and 
Taylor,  1995;  and  Jih’s  talk  at  this  SSA  for  some  aspects  of  this  issue.)  This  is  a  typical  example  of 
uncertainties  in  our  physical  understanding  of  regional  phase  excitation  due  to  the  (current)  inseparability 
of  source  and  path  effects. 

A  basic  research  program  should  also  include  work  on  a  predictive  capability  that  can  model  fully 
3-D  elastic  propagation  efficiently  in  the  high-frequency  (1-10  Hz)  regime  out  to  regional  distances  (1000 
km  and  further).  The  study  of  the  blockage  of  other  regional  phases  besides  Lg  is  also  important,  espe¬ 
cially  as  it  relates  to  location  and  discrimination. 

We  think  an  applied  program  should  complement  the  basic  research  studies  with  projects  aimed  at 
empiricaly  mapping  blockage  in  regions  of  CTBT  interest.  The  could  include  projects  using  both  histori¬ 
cal  (such  as  the  ILPA  dataset)  and  current  data  sets.  This  type  of  work  will  be  important  in  completing 
the  calibration  of  different  regions  using  regional  seismic  data.  In  addition,  the  applied  program  is  the 
appropriate  place  to  conduct  benchmark  tests  of  different  numerical  modeling  techniques. 

Both  the  basic  and  applied  research  programs  could  involve  field  studies;  however,  the  design  and 
purpose  of  these  experiments  should  be  very  carefully  planned.  Unless  the  proposed  experiment  is 
carefully  designed  to  solve  some  specific  scientific  puzzle,  it  might  end  up  to  be  a  nonproductive  training 
experiment.  During  1995  DoE  conducted  five  field  experiments  (Stump  et  al.,  1996;  DoE,  1996),  and 
the  Air  Force  (PL  and  AFOSR)  also  funded  at  least  two  experiments,  in  our  opinion  it  may  be  preferable 
to  digest  the  data  from  these  experiments  somewhat  before  more  field  experiments  are  conducted. 
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Possible  basic  and  applied  types  of  field  experiments: 

1.  A  6.1  experiment  should  be  designed  to  solve  some  fundamental  scientific  puzzle  or  issue  that 
is  relevant  to  Air  Force’s  treaty  monitoring  need.  Currently  the  most  important  outstanding  question  is 
the  relationship  between  near-source  media  and  the  mechanics  of  sources  to  radiated  seismic  signals. 
This  type  of  question  cannot  be  answered  by  refraction  surveys,  simply  deploying  more  sensors  or  mak¬ 
ing  more  field  trips. 

A  good  example  of  a  relevant  6.1  experiment  is  the  DoE  field  program  led  by  Brian  Stump’s  group 
at  Los  Alamos  National  Laboratory.  During  1995,  Stump  et  al.  (1996)  conducted  a  series  of  field  exper¬ 
iments  with  very  specific  goals  attached  to  each  experiment.  For  instance,  his  “Source  Geometry’’ 
experiment  was  designed  to  compare  spherical  (typical  of  tamped  underground  nuclear  explosion)  and 
cylindrical  explosions  (typical  of  mining  explosions)  and  to  develop  monitoring  techniques  for  both  types. 
In  the  “Black  Thunder”  experiment.  Stump’s  research  goal  was  to  document  the  effect  of  blasting  prac¬ 
tice  on  observable  seismic  signals  (see  also  DoE,  1996). 

Although  a  basic  research  program  like  AFOSR’s  is  the  ideal  vehicle  to  fund  experiments  of  the 

6.1  type,  such  experiments  (like  Stump’s)  are  prohibitively  expensive.  As  a  result,  only  very  few  Pis  can 
be  funded  to  conduct  this  type  of  research.  What  the  6.1  community  should  do  is  perhaps  to  team  with 
Stump  and  analyze  his  data  jointly,  rather  than  trying  to  duplicate  his  experiments  at  the  expense  of 
other  fundamental  research.  Perhaps  the  funding  agencies  should  make  a  major  field  experiment  every 
two  or  three  years  a  central  element  of  all  of  the  agencies’  programs.  This  would  take  major  coordina¬ 
tion  efforts  among  AFTAC,  DARPA,  DoE,  DNA,  Phillips  Laboratory,  and  AFOSR,  plus  potential  partici¬ 
pants  from  foreign  countries,  but  the  experiment  could  sen/e  to  develop  some  commonality  among  the 
investigators. 

2.  A  6.2/6.3  type  of  field  program  should  include  reconnaissance  surveys  for  regions  of  monitoring 
concern  that  have  not  been  explored  well  before.  Agencies  that  manage  6.2/6.3  programs  should  fund 
studies  that  attempt  to  gather  useful  data  to  form  an  usable  archive.  For  this  type  of  experiment,  the 
selection  of  site  is  critical,  and  the  basic  research  program  should  therefore  not  be  obligated  to  fund  this 
type  of  experiment  In  fact,  if  a  PI  is  interested  in  “calibrating”  a  region  of  high  monitoring  concern 
through  a  series  of  field  experiments,  the  relevant  agencies  who  would  use  the  results  had  better  be 
consulted  in  advance. 

The  worst  scenario  is  for  a  PI  to  be  funded  under  6.1  to  conduct  an  experiment  with  no  scientific 
goal  in  mind  or  funded  under  6.2  to  conduct  experiments  in  an  area  of  no  6.2  relevance  at  all.  Such 
experiments  are  neither  a  6.1  nor  6.2  project  -  they  serve  at  best  as  a  training  trip  for  the  technicians. 
Such  experiments  bear  no  relevance  to  the  Air  Force’s  mission  and  should  not  be  funded  by  the  6.1  or 

6.2  programs,  especially  when  the  budgets  of  these  programs  are  already  extremely  tight. 
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Lg  Propagation-Path  Barriers  in  the  Eurasian  Continental  Craton  -  Possible  Shallow  Crust 
Explanations 

Baumgardt,  D.  R.,  ENSCO,  Inc.,  5400  Port  Royal  Rd.,  Springfield,  VA  22151,  doug@ensco.com. 

Investigations  of  regional  seismic  phase  propagation  in  Eurasia,  including  Scandinavia,  Europe, 
Russia,  and  China,  have  revealed  numerous  instances  of  anomalous  attenuation,  or  "blockage ",  of  the 
regional  shear  phase,  Lg.  Because  many  regional  seismic  discriminants  used  to  monitor  comprehensive 
test  ban  treaties  rely  on  Lg  amplitude  and  spectral  measurements,  it  is  important  to  understand  the  origin 
of  sudden  Lg  attenuation. 

Lg  barriers  in  the  Eurasian  continental  cratons  are  related  to  laterally  heterogeneous  shallow  cru¬ 
stal  structure.  Such  heterogeneities  appear  to  distort  the  crustal  waveguide  for  Lg  and  thus  produce  a 
sudden  blockage  of  Lg  propagation  through  the  anomalous  region.  In  the  continents,  Lg  blockages  have 
been  observed  for  deep  sedimentary  basins  which  act  as  Lg  "sinks",  i.e. ,  they  trap  Lg  waves  which  pass 
through  due  to  strong  impedance  contrasts  between  the  low-velocity  sediments  in  the  basins  and 
granitic-velocity  rocks  outside  of  the  basins.  Crustal  thickness  variations  (e.g.,  crustal  pinchouts)  do  not 
seem  to  be  implicated  in  these  blockages  since  the  crustal  thickness  does  not  vary  by  more  than  10  km 
in  these  regions,  even  under  mountainous  regions  like  the  Urals.  Examples  of  such  regions  include  the 
Barents  Sea,  Pechora,  Caspian  Basins,  and  sedimentary  basins  around  the  Urals  in  Russia. 

in  China,  strong  blockage  of  high-frequency  (>1  Hz)  Lg  has  been  observed  for  propagation  paths 
crossing  southern  Tibet,  although  Lg  does  appear  to  propagate  efficiently  across  northern  Tibet.  More¬ 
over,  low-frequency  (<1  Hz)  Lg  does  not  seem  to  be  blocked  in  any  of  these  regions.  Lg  blockage  in 
Tibet  has  usually  been  attributed  to  anomalous  crustal  thickness  and  high  anelastic  attenuation. 
Although  Tibet  is  ringed  by  sedimentary  basins,  such  as  the  Tarim  Basin  in  the  north,  the  sediment 
depths  are  shallow  and  do  not  seem  to  strongly  block  Lg  in  the  north.  However,  since  many  of  the 
blocked  Lg  paths  in  southern  Tibet  cross  the  Tsangpo/Indus  suture  zone,  a  region  of  anomalous  geologi¬ 
cal  structure,  it  is  possible  that  this  suture  zone  may  be  a  thin  boundary  in  southern  Tibet  that  blocks 
high-frequency  but  not  low-frequency  Lg  waves. 

These  results  suggest  that  lateral  heterogeneity  in  shallow  crustal  structure  may  be  a  major  cause 
of  high-frequency  Lg  blockage  in  the  continents.  More  observations  from  new  stations  now  appearing  in 
Eurasia  will  be  needed  to  more  thoroughly  map  out  regions  of  anomalous  Lg  block^e  in  order  to  explain 
their  origin  and  to  properly  calibrate  regional  discriminants  which  utilize  measurements  of  the  Lg  phase. 


Anomalous  Surfoce  Wave  Propagation  across  the  South  Caspian  Basin  and  the  Block¬ 
age  of  regional  Seismic  Phases 

K.  Priestley,  Bullard  Laboratories,  Univ.  of  Cambridge,  Cambridge  CBS  OEZ,  United  Kingdom, 
keith@esc.cam.ac.uk,  C.  Shultz,  and  H.  Patton,  Division  of  Earth  Sciences,  Lawrence  Livermore 
national  Laboratory,  CA  94550. 

The  crust  and  upper  mantle  structure  of  the  south  Caspian  Basin  is  enigmatic.  Early  Soviet  studies 
show  that  the  crust  of  the  basin  consists  of  two  layers:  a  thick  sedimentary  section  (15-25  km)  with  low 
P-wave  velocity  (3.5-4.0  km/s)  overlying  a  12-18  km  thick  basaltic  lower  crust  The  study  of  Kadinsl^- 
Cade  et  al.  (1981)  demonstrated  that  the  seismic  phase  Lg  is  largely  blocked  for  paths  crossing  the 
south  Caspian  Basin.  New  seismic  data  shows  that  the  south  Caspian  Basin  also  severely  disrupts  low 
frequency  (0.017  -0.10  hiz)  fundamental  mode  surface  wave  trains.  The  effect  is  observed  for  surface 
waves  propagating  along  both  east-to-west  and  west-to-east  great  drde  paths  showing  that  this  is  not  a 
site  or  instrumental  effect.  We  model  the  response  of  the  surface  wave  to  this  low-velocity  sediment, 
deep  basin  structure  and  crustal  thinning  using  a  hybrid  locked-mode/finite-difference  approach.  We 
demonstrate  that  much  of  the  observed  surface  wave  train  degradation  can  be  modeled  with  the  2-D 
basin  structure.  Finally  we  demonstrate  how  a  simple  model  of  the  structure  in  the  Caspian  region  can 
help  identify  blockage  mechanism. 
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Numerical  Modeling  of  Rg  Blockage  and  Comparison  of  Bg-to-Lg  Conversion  Mechanisms 
due  to  Scattering  and  Incomplete  Dissipation 

Jih,  R.-S.,  PUGPE,  Hanscom  AFB,  MA  01731,  jih@snappy.plh.af.mil;  McLaughlin,  K  L,  S-CUBED,  La 

Jolla,  CA;  Dainty,  A.  M.,  PL/GPE,  Hanscom  AFB,  MA;  and  Harkrider,  D.  G.,  Boston  College  @ 

PUGPE,  Hanscom  AFB,  MA. 

The  presence  of  a  short-period  fundamental-mode  Rayleigh  wave,  Rg,  is  an  excellent  diagnostic  for 
shallow  source  depth.  However,  Rg  is  susceptible  to  non-uniformities  of  the  near-surface  waveguide, 
and  hence  absence  of  Rg  does  not  imply  a  deep  source.  Without  a  predictive  capability,  the  us^e  of 
Rg  will  be  limited  to  well-caEbrated  source  regions  and  propagation  paths.  Study  of  Rg  propagation 
(scattering,  attenuation,  etc.)  would  naturally  help  constrain  source  depth  and  thereby  help  event 
identification  also.  In  this  study,  intensive  linear  finite-difterence  [LFD]  calculations  are  conducted  to 
explore  the  physical  mechanisms  of  Rg  blockage  with  emphasis  placed  on  quantification  of  the  scattering 
effects  due  to  rough  topography  and  random  heterogeneity.  In  addition,  a  new  Rg  -to-  Lg  conversion 
mechanism  due  to  incomplete  dissipation  by  shallow,  strong  attenuative  layers  is  presented  and  the  rela¬ 
tive  effectiveness  of  all  three  Rg-related  mechanisms  are  compared.  For  models  embedded  with  shallow 
random  heterogeneity,  the  RMS  velocity  fluctuation  correlates  very  well  with  the  Rg  transmission 
coefficient.  For  1  Hz  Rg,  a  2%-5%  variation  in  the  velocity  leads  to  an  equivalent  spatial  Q  value  of 
several  hundreds  or  larger,  regardless  of  which  of  the  three  commonly  used  random  media  is  embed¬ 
ded.  Rough  topography  typically  results  in  a  Q  value  ranging  from  10  to  100,  which  is  approximately 
equivalent  to  a  random  medium  with  velocity  variation  larger  than  10%.  Incomplete  dissipation  of  Rg 
waves  produces  a  very  simple  wave  field,  with  almost  all  of  the  undissipated  energy  propagating  laterally 
towards  the  forward,  and  hence  postcritical,  directions.  On  the  other  hand,  the  scattering  by  shallow 
heterogeneity  or  rough  topography  generates  a  rather  complicated  wave  field,  with  a  significant  fraction 
of  the  scattered  energy  going  steeply  downward.  In  terms  of  energy  partitioning  between  converted  P 
and  S  waves,  the  incomplete  dissipation  overwhelmingly  directs  more  energy  into  S  (Lg)  waves,  as  com¬ 
pared  to  the  remaining  two  mechanisms.  For  whichever  of  the  three  Rg-bloddng  mechanisms  is 
invoked,  the  S  waves  converted  from  Rg  always  dominate  the  whole  scattered  field. 
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Pn  amplitude  variation  along  the  northern  foot  of  the  Tienshan  mountains 

Xie,  J.,  Lamont-Doherty  Earth  Observatory,  Columbia  University,  Route  9W,  Palisades,  NY  10964;  914- 

365-8553;  xie@lamont.ldgo.columbia.edu. 

The  term  "blockage"  of  regional  waveforms  has  been  used  to  describe  drastic  amplitude  reduc¬ 
tions  of  regional  waveforms.  This  term,  however,  is  typically  used  in  a  qualitative  and  relative  sense. 
For  example,  blockages  of  Lg  or  Sn  phases  require  less  amplitude  reduction  than  that  of  the  primary  (P) 
phase  since  Lg  or  Sn  are  preceded  by  coda  of  other  phases  which  may  have  considerable  amplitude,  in 
reality,  the  amplitude  of  P  waves  often  vary  much  more  drastically  than  that  of  the  Lg  phase,  and  in  a 
much  more  unpredictable  manner.  We  present  a  case  in  which  the  Pn  phase  from  Lop  Nor  explosions 
undeigoes  a  drastic  variation  in  amplitude  (by  a  factor  of  30)  aaoss  the  Kyighizstan  Network  (KNET). 
This  network  is  deployed  along  the  northern  foot  of  the  Tienshan  mountains  and  has  a  relatively  small  (« 
100  km)  aperture.  The  amplitude  variation  primarily  occurs  at  lower  (s  1  Hz)  frequencies,  and  is  accom¬ 
panied  by  anomalous  polarizations.  At  station  KZA  where  the  largest  Pn  amplitude  is  observed,  Pn  is 
polarized  to  N60°W,  about  30°  away  from  the  great  circle.  The  cause  of  the  amplitude  and  polarization 
anomalies  are  likely  to  be  some  deep  3D  stnictural  variations,  such  as  Moho  topography. 

Unlike  the  Pn  amplitude,  the  Lg  amplitude  is  quite  stable  aaoss  the  KNET.  This  suggests  that  the 
amplitude  of  Pn  is  much  more  sensitive  to  smaller  scale  variations  of  crustal  structure  than  that  of  the  Lg 
phase. 
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Regional  variations  In  Lg  observed  and  synthesized  at  the  CNET  and  KNET  arrays 

Anderson,  T.  S.,  and  V.  F.  Cormier,  Dept,  of  Geology  and  Geophysics,  University  of  Connecticut,  Storrs, 

CT  06269,  tom@geol.uconn.edu,  cormiet@geol.uconn.edu. 

Lg  coda  seismograms  are  synthesized  as  multiple  SmS  waves  in  3-D  models  of  the  crust  and 
upper  mantle  by  dynamic  ray  tracing  and  summation  of  Gaussian  beams.  Results  are  compared  with 
data  from  the  CNET  (Caucasus)  and  KNET  (Kyrgyzstan)  arrays,  located  in  regions  having  significant 
crustal  variations  along  the  borders  of  the  former  Soviet  Union.  Ray  diagrams  predict  weak  regional 
variations  in  Lg  efficiency  associated  with  weak  gradients  in  Moho  topography  in  the  Caucasus  region.  In 
this  region,  synthetic  seismograms  predict  crustal  thickness  variations  will  have  little  effect  on  variations 
seen  across  the  network.  Thus,  effects  other  than  crustal  thickness  variations  must  be  sought  to  explain 
strong  regional  variations  in  Lg  waveforms  observed  from  different  azimuths  at  CNET.  Ray  diagrams 
predict  much  stronger  regional  variations  in  Lg  efficiency  in  the  Kyrgyzstan  region.  These  variations 
correlate  well  with  strong  gradients  in  Moho  topography  along  the  Hindu  Kush  and  Pamir  mountain 
ranges.  Detailed  variations  of  waveforms  across  the  arrays  are  consistent  with  Lg  efficiency  being  pro¬ 
portional  to  the  length  and  number  of  times  SmS  ray  paths  traverse  regions  of  strong  Moho  gradient. 
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Regional  Phase  Propagation  in  Northern  Africa  and  the  Mediterranean 

McNamara,  D.  E.  W.  R.  Walter,  C.  Schultz  and  P.  Goldstein,  Livermore  National  Laboratory,  Earth  Sci¬ 
ences  Division,  L-205,  Livermore,  CA  94551. 

The  Mediterranean  region  is  characterized  by  abundant  seismicity  in  a  complex  tectonic  environ¬ 
ment  resulting  from  the  convergence  of  the  African  and  Aslan  plates.  The  purpose  of  this  study  is  to 
analyze  the  propagation  characteristics  of  the  high  frequency  regional  phases  Pn,  Pg,  Sn  and  Lg  within 
North  Africa  and  across  the  Mediterranean  Sea.  We  have  collected  regional  seismograms  from  North 
African  earthquakes  recorded  at  MEDNET,  GEOSCOPE  and  IRIS  stations  in  the  area.  A  first  look  at  the 
data  indicates  that  Lg  propagates  efficiently  to  distances  of  2000  km  or  more  within  Northern  Africa  but 
is  effectively  blocked  for  most  paths  crossing  the  Mediterranean.  The  picture  for  Sn  is  more  complicated 
and  shows  no  dear  regional  pattern.  However  unlike  Lg,  Sn  is  observed  on  some  paths  crossing  the 
Mediterranean  Sea.  These  S-wave  phases  (Sn  and  Lg)  are  particulariy  important  for  algorithms  that 
discriminate  small  magnitude  (Ms  <  3.5)  underground  explosions  from  earthquakes.  Many  studies  have 
shown  that  changes  in  crustal  thickness  and/or  upper  mantle  velodty  can  change  or  even  block  the 
transmission  of  Sn  and  Lg,  greatly  complicating  the  problem  of  identifying  dandestine  explosions  in  such 
regions.  Guided  by  our  empirical  data,  we  make  use  of  2-D  modeling  codes  to  better  quantify  the  types 
of  structures  and  lengths  of  paths  within  those  structures,  that  block  these  regional  phases. 

This  research  was  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  Lawrence  Livermore 
National  Laboratory  under  contract  number  W-7405-ENG-48. 


Interpretation  of  Crustal  Phase  Characteristics  in  Iran  and  the  Surrounding  Region 
Determined  from  ILPA  Data 

Sweeney,  J.  J.,  Treaty  Verification  Program,  Lawrence  Livermore  National  Laboratory,  P.O.  Box  808,  L- 
208,  Livermore,  CA,  94550,  sweeney3@llnl.gov 

The  presence  or  absence  of  crustal  phases  such  as  Sn  and  L^  at  the  Iranian  Long  Period  Array 
(ILPA)  show  a  complex  pattern  in  Iran  and  the  surrounding  region.  A  few  degrees  of  azimuth  may  mark 
the  boundary  between  an  area  where  a  phase  is  present  and  where  it  is  absent.  Previous  work  has 
identified  the  southern  Caspian  region  as  an  area  of  poor  Sn  propagation.  This  study  confirms  that  result 
and  also  identifies  southeastern  Makran,  part  of  the  Iran  plateau,  and  part  of  the  Z^ros  region  as  areas 
where  the  Sn  phase  is  absent.  In  general,  lack  of  an  Sn  phase  can  be  associated  with  paths  aossing 
regions  of  Ckratemary  and  Tertiary  volcanism.  Lg  propagation  is  present  over  most  of  Iran  with  the 
exception  of  the  Zagros  region  between  the  Oman  line  and  Qatar.  This  is  an  area  with  extensive  mobile 
infra-Cambrian  salt.  Other  regions  identified  with  no  Lg  phase  are  the  Caspian  and  Black  Sea  areas, 
which  have  been  noted  by  other  workers. 

This  work  was  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  Lawrence  Livermore  National 
Laboratory  under  contract  No.  W-7405-ENG-48.  Multimax  Corp.  supplied  the  data  for  this  study. 
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